Fluctuation theorem for entropy production during effusion of a relativistic ideal gas 
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The probability distribution of the entropy production for the effusion of a relativistic ideal gas 
is calculated explicitly. This result is then extended to include particle and anti-particle pair pro- 
duction and annihilation. In both cases, the fluctuation theorem is verified. 
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I. INTRODUCTION 



The fluctuation theorem [l|, [3, H, 01 can be seen as a 
natural extension to far from equilibrium situations of the 
pioneering work of Onsager on the reciprocity relations. 
The theorem has been proven in a variety of settings, in- 
cluding thermostated systems [1, H, 01 , Markov processes 

111,111, 



17| and 



eluding thermosi 

d, [^, Uo, EH , classical Hamiltonian dynamics 
quantum dynamics, quantum field theory (l5l.[l^ 
relativity [ISi] . The wide range of validity of the theorem 
reflects the fact that it derives from two basic ingredients 
of the fundamental laws of nature: time reversibility and 
time translational invariancc. The effusion of gases pro- 
vides a simple setting in which the fluctuation theorem 
can be verified on the basis of freshmen's physics [l^, |23] ■ 
The effusion of radiation (photons) can be discussed on 
a similar basis [111]. In this paper, we show that the fluc- 
tuation theorem is also valid for the effusion of classical 
relativistic particles, including particle and anti-particle 
pair production and annihilation. 



II. FLUCTUATION THEOREM FOR EFFUSION 
OF A RELATIVISTIC IDEAL GAS 



Two relativistic ideal gases, A and B, are initially in 
equilibrium at given densities pi and temperature Ti, 
i S {A, B} while separated by a common adiabatic wall. 
At time t = 0, a small hole with surface area a is opened, 
allowing an exchange of particles and energy between the 
two gases during a (fixed) time interval of length t. We 
assume the presence of an energy filter, so that only par- 
ticles with total energy in the range [Ei, E2] can cross. 
The dimensions of the hole are supposed to be small com- 
pared to the mean free path of the gas particles, and the 
reservoirs sufficiently large so that the equilibrium state 
in both parts is not affected by the exchange. In par- 
ticular, the temperatures and densities in both reservoirs 
remain constant. Hence, the change in entropy in the 
total system, upon a total transfer of energy AC/ and of 



particles AiV from A to B during is given by: 
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We introduced, in accordance with the definitions from 
irreversible thermodynamics, the following thermody- 
namic forces (affinities) for energy and particle flow, re- 
spectively: 

- 7f , - 7f 7T^- i^) 

A derivation of the chemical potential /ii of a relativis- 
tic ideal gas is given in Appendix \^ and leads to the 
following explicit expression for An'- 

' PaTb K2 {mc^ /kTB)' 



An — k log 
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where K2 is the (second) modifled Bessel function of the 
second kind. In writing down these results, we use the 
convention that the internal energy U of the gas com- 
prises both the kinetic energy and the rest energy of the 
particles {mc^). Note that the thermodynamic force An 
diverges when one of the (infinitely large!) reservoirs is 
empty, so that free effusion into unlimited space, imply- 
ing an infinitely large entropy production, corresponds to 
a singular limit. 

We will show below that the joint probability density 
Pt (AC/, AN) satisfies the following detailed fluctuation 
theorem: 



Pt{AU,AN) 
Pt{-AU,^AN) 
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with AS' given by Eq. ([T]). This encompasses as a special 
case the more familiar form of the fluctuation theorem: 



Pt{AS) 



Pti-AS) ■ 

Since the increments of AU and AA^ in successive time 
intervals are independent, one has that: 

J2 /"e-(^^^^+^"^^)Pt(AC/,AiV)dAC/ = 

AN '' 

g-tQ(A[/,Ajv)_ 
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where Q{Xu,^n) is the cumulant generating function. 
In terms of the latter, the detailed fluctuation theorem 
reads: 



Q{\u,Xn) = QiAu/k- \u,AN/k- Xn)- (7) 



III. MASTER EQUATION AND CUMULANT 
GENERATING FUNCTION 



During a small time interval dt, the contributions to 
quantities AU and AN result from individual particle 
transport across the hole. Following basic arguments 
from kinetic theory of gases the corresponding probabil- 
ities per unit time, Ta^b{E) and Tb^a{E), to observe 
a particle with energy E crossing the hole from i ^ j \s 
given by (cf. Appendix [B]) : 



api {mc'^/kTi)e '■^i 
4c2(mc)3 K2{mc^/kTi) 
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The probability density Pt{AU, AN) thus obeys the fol- 
lowing Master equation: 
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PtiAU, AN) [ ' [Ta^b (E) + Tb^a (E)] dE, (9) 
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subject to the initial condition Po{AU, AN) — 
S{AU)San,o- 

The cumulant generating function Q{Xu, Aat) is obtained 
by multiplying the master equation with exp{— A£/AJ7 — 
XjsfAN} and subsequent integration over AU and sum- 
mation over AN respectively. This leads to the following 
analytical result: 



Q{Xu, Xn) = IaIO.O) - IaIXu, Xn) 

+ Ib{0,0)-Ib{-Xu,-Xn), (10) 

where we define: 



liiXu, Xn) = 
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The evaluation of the integral yields: 
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Eq. (jlOp. together with Eq. (fT^ . is the central result of 
this paper. A tedious calculation, using the explicit ex- 
pressions for the thermodynamic forces given in Eq. ([3]), 
reveals that this expression indeed verifies the detailed 
fluctuation theorem Eq. ([7]). 

The analysis becomes somewhat simpler in the extreme 
relativistic limit mc^/kT ^ 1. In this case the rela- 
tivistic energy becomes H{p) = cp and the results from 
Appendices [X] and |B] simplify: 
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The same expressions result of course also from (jASp . 

and ([3]) respectively by using Ei ^ mc^ and 
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K2{u) ~ 2/^2 for u < 1. Eq. ^ simplifies to 



Ii{Xu,XN) 



apiC 



4(1 + XukTi)^ 
and the cumulant generating function becomes 
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The fluctuation theorem ([7]) is easily verifled. 



IV. CUMULANTS 

A Taylor expansion of the cumulant generating func- 
tion gives immediately the expressions of the joint cumu- 
lants Kij (of power i in the energy flux and power j in 
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the particle flux): 

q{xu,xn) = (19) 

As the higher order expressions quickly become unwieldy, 
we will only present the first order results, namely kiq = 

I 



(At/) and kqi = (AN). These results will be used in the 
next section to validate the Onsager symmetry. 
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V. ONSAGER RELATION 

The average entropy production during the effusion 
process is obtained by averaging Eq. ([1]) and taking the 
time derivative: 



dt 



(AS) ^AuJu + AnJn, 



(22) 



where we introduce the macroscopic fluxes jTiy and JTat 
corresponding to energy and particle transport: 



Ju = j^{AU) 



J^^-{AN). (23) 



These fluxes in general depend nonlinearly on the ther- 
modynamic forces. Close to equilibrium, when the ther- 



modynamic forces are small (that is, for small tempera- 
ture and density difference), the expressions can be lin- 
earised, and written in compact matrix notation 



(24) 



with J = {Ju,JnV, a = {Au,AnV and L the On- 
sager matrix : 



L = 



(25) 



The following explicit expression for these Onsager coef- 
ficients are obtained from Eqs. (|20p and pT|) : 
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The symmetry property of the matrix is a direct conse- 
quence of the symmetry relation Eq. ([7]) of the cumulant 
generating function. 

These expressions significantly simplify when considering 
an infinitesimal energy window (dE <^ Eq): 



El = Eq- dE/2 ; 
E2 = EQ + dE/2 . 



(27) 



In this case, there is a strong coupling between the parti- 
cle and energy current: every particle that crosses carries 
a energy Eq and so J^u — EqJ^. This is reflected in the 
Onsager matrix, after expanding the to first order in 
dE: 
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The determinant of the Onsager matrix is identically 
zero, as both rows are linear dependent. This opens up 
the possiblity, in principle, to use the system as a ther- 
modynamic engine attaining Carnot or Curzon-Ahlborn 
efficiencies (22, \2A . 



VI. PAIR PRODUCTION 

One important aspect of relativistic gases, which we 
have not considered so far, is the production and anni- 
hilation of particle-antiparticle pairs ,24]. This must be 
included in a consistent treatment of relativistic effusion 
by considering the corresponding exchange of antiparti- 
cles between the two compartments. 
In our framework, in which both reservoirs independently 
are (and remain) in equilibrium, this implies the exis- 
tence of a second ideal gas in each reservoir consisting of 
anti-particles. Because of the production and annihila- 
tion, both particles and anti-particles will be created and 



destroyed. This process can be represented by a chemical 
reaction of the following form: 



(29) 



(x represents a particle, x the corresponding anti- 
particle). Equilibrium is achieved when the two chemical 
potentials, p, and /I, satisfy: 



/X + /i = 0. 



(30) 



For a given temperature this relation determines the den- 
sity p of the anti-particles. Using Eq. (jA4[) gives: 
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The same result can be obtained by a direct calculation 
of the grand canonical partition function of the two ideal 
gases, taking into account that for pair production and 
annihilation the difference between the number of parti- 
cles and antiparticles is conserved [25j . 
Note that fluctuations of (anti-)particles number and en- 
ergy around their equilibrium values can be neglected, as 
both reservoirs are assumed to be infinitely large. 
The exchange of antiparticles through the opening be- 
tween volumes A and B will contribute to the entropy 
change AS". Since the chemical potential of the antipar- 
ticles is just the negative of the chemical potential of the 
particles we find instead of ([1]) 



AS = AuAU + AnAN + AnAN 
= AuAU + ANiAN - AN), 



(32) 



where Au and An{— ^Aj^) are given by ^ and ^ 
respectively. To account for the fluctuations in AS" due 
to the exchange of antiparticles the joint probability dis- 
tribution has to be augmented to P{AU, AN, AN), and 
the corresponding cumulant generating function is given 
by: 



-i\uAU+\NAN+\g,AN)\ 
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The detailed fluctuation theorem acquires the form The time evolution of P{AU, AN, AN) is now governed 

by the master equation: 

Q{Xu, Aat, A^) = 

Q{Au/k~ \u,AN/k~ XN,AN/k- Xj^). (34) 



^Pt{AU,AN,AN)^ f ' [TA^B{E)Pt{AU - E,AN -l,AN) + fA^BiE)Pt{AU - E,AN,AN -1) 

+ TB^A{E)PtiAU + E,AN +l,AN) +TB^A{E)PtiAU + E,AN, AN + I) 
-Pt{AU, AN, AN) {Ta^b{E) + fA^B{E) + Tb^a{E) + Tb^a(S))] dE, (35) 



where the rates Ti-,j{E) and Ti^j{E) are defined by 
(|15|1 with Pi replaced by pi in the latter case. Similar 
to Eq. (fTU]) we now find: 



Q{Xu,Xn,Xj^) = Ia{0,0)- Ia{Xu,Xn)+Ia{0,0) 
- Ia{Xu,Xn) + Ib{0,0) - Ib{-Xu,-Xn) 

+ /b(0,0)-/b(-Ac/,-A^) , (36) 

with Ii{Xu, Xn) given by Eq. For li the density is pi. 
Using these expressions in (j36|) it is again straightforward 
to verify ([M)) . 

The inclusion of pair production thus opens an ad- 
ditional channel for entropy production. However, due 
to the relations Eqs. I|30ll3ip between the chemical po- 
tentials of particles and antiparticles and their respec- 
tive concentrations in equilibrium the fluctuation theo- 
rem nevertheless remains valid. 

In the above, we have considered only one kind of 
particle-antiparticle pair. The extension to different 
kinds of pairs is straightforward. The contribution of 
the photons produced in the annihilation processes to 
the entropy production may be included along the lines 
of El. 



If the gas is enclosed in a volume V, the expression for 
Z becomes: 
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where K2 is the second modified Bessel function of the 
second kind. 

Given Z, it is now standard procedure to obtain the ther- 
modynamic quantities of interest. The chemical poten- 
tial p can be derived from the Helmholtz free energy F 
(p^N/V): 



M = 



OF _ rATrimc)^ K2imc^/kT)\ 
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with F = —kTlnZ given by (in the thermodynamic 
limit): 
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APPENDIX B: TRANSITION RATES 



APPENDIX A: THERMODYNAMICAL 
PROPERTIES OF A RELATIVISTIC IDEAL GAS 



The basic thermodynamic properties follow directly 
from the partition function Z , which for a gas of N iden- 
tical non-interacting particles is given by: 
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The Hamiltonian of a free relativistic particle is given by 
the usual expression of the relativistic energy: 



(A2) 



The calculation of the transition rates as given in 
Eq. ([8]), is similar to that of the pressure in textbooks 
on kinetic theory. Ti^j{E)dEdt is the probability to ob- 
serve a particle with energy (kinetic -I- rest energy) in the 
range [E, E + dE] , crossing the hole from i to j during 
a time interval dt. The corresponding momentum of the 
particle must be in the range [p,p + dp] with: 



P 



m^c^ and dp 
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Since the gas is in equilibrium, its momentum distribu- 
tion has the usual Boltzmann form: 



^.■K(rncf'K2(jn(? jkTi) 



-H{p)/kT, 
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with H(v) given in Eg. ( jA2[) . As the particle (with veloc- Finally, replacing v = pc/ y^p^ + mc^ and integrating 
ity v) must be able to reach the hole during the time inter- yields the expression given in Eq. ([5]). 
val dt, it must be located inside the volume av cos 9dt {9 is 
the usual spherical coordinate, taking the z-axis perpen- 
dicular to the common wall between the two reservoirs 
and pointing from i to j). Since the particle is moving 
towards the hole we have < 9 < 7r/2 . Adding all con- 
tributions from the different directions 9 and tp leads to 
the following result: 

T^^j{E)dEdt 

= / p^ siT\9dpd9dLp Pi4>i{p) av COS 9dt. (B3) 
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